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We study the previously conjectured form of the physical electron propagator and its allegedly 
Luttinger type of behavior in the theory of the pseudogap phase of high-temperature copper-oxide su- 
perconductors and other effective QED-like models. We demonstrate that, among a whole family of 
seemingly gauge-invariant functions, the conjectured "stringy ansatz" for the electron propagator is 
the only one that is truly invariant. However, contrary to the results of the earlier works, it appears to 
have a negative anomalous dimension, which makes it a rather poor candidate to the role of the phys- 
ical electron propagator. Instead, we argue that the latter may, in fact, feature a "super-Luttinger" 
behavior characterized by a faster than any power-law decay Gphys{x) oc exp(— const In^ |a;|). 
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O ■ I. INTRODUCTION 

(N : 

. In a generic (1 + l)-dimensional (hereafter 2D) many-fermion system, an arbitrarily weak short-range repulsive 
interaction is known to completely destroy the conventional Fermi liquid, thereby giving rise to the so-called Luttinger 
: behavior. As one of the hallmarks of the Luttinger regime, the electron propagator exhibits an algebraic decay with 
' distance which is controlled by a non-universal exponent and appears to be faster than in the non-interacting case. 

In higher dimensions (3D and 4D) the Fermi liquid is believed to be more robust, although it is not expected 
' to remain absolutely stable. In the absence of a clear-cut evidence, however, the possibility of a non-Fermi liquid 
, behavior in D > 2 has recently become the subject of an intense debate. 

' While in the case of short-ranged repulsive interactions any departures from the Fermi liquid are likely to be limited 
to the infinitely strong coupling regime, the long-ranged forces are considered to be capable of destroying the Fermi 
\l liquid even at finite coupling. One of the most extensively studied such examples, the finite density system of non- 
relativistic massive fermions which are minimally coupled to an abelian gauge field, was indeed shown to manifest a 
distinctly non- Fermi liquid behavior, although the latter turned out to be rather different from the Luttinger one [1]. 
More recently, there has been a renewed interest in the 3D relativistic counterpart of this model governed by the 
' conventional QED3 action 
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^ : for the A^-flavored Dirac fermions which are described by a reducible (four-component) representation corresponding 
' to the choice of the 7- matrices 7^^ = ct^ 173 constructed from the triplet of Pauli matrices. The use of the 
• ^ : parity-even representation guarantees that radiation corrections do not generate any parity-breaking Chern-Simons 
^\ terms. 

^ I Among the previously discussed examples of the 3D condensed matter systems that support the Dirac-like low- 
" " ' energy excitations and allow for such an effective description are the so-called flux phase in doped Mott insulators 
[2,3] and d-wave superconductors with strong phase fluctuations [4-7] proposed as an explanation of the pseudogap 
and insulating (spin and/or charge density wave) phases of the high-Tc cuprates. Also, the non-Lorentz- invariant 
version of QED3 was shown to provide a convenient description of the normal semimetalic state of highly oriented 
pyrolytic graphite [8-10]. 

The number A^ of the fermion flavors depends on the problem in question, although it is not necessarily equal 
to the number of different conical Dirac points in the bare electron spectrum. In the abovcmcntioncd condensed 
matter-inspired QE D-likc models [2-10], A^ = 2 is a number of the electron spin components, while the number of 
the conical points turns out to be either four [4-7] or two [8-10]. 

In some of the strongly correlated electron systems described by the effective QED-like theories, such as the problem 
of graphite, the theory is formulated in terms of the electrons themselves which enables one to directly compute various 
experimentally relevant observables [9,10]. 

However, in other systems the effective gauge flelds serve as a convenient representation of such bosonic collective 
modes as spin or pairing fluctuations, while the Dirac fermions correspond to some auxiliary fermionic excitations, 
such as spinous [2-4], "topological" fermions [5-7], and so forth. 

In the latter case, a generic quantum mechanical amplitude written in terms of the Lagrangian fermions turns out 
to be gauge-dependent, while all the physical observables, which experimental probes can only couple to, must be 
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manifestly gauge-invariant. In light of that, it is absolutely imperative to establish the correct form of the physical 
electron propagator in terms of the Lagrangian fermions. Until this task is accomplished, no theoretical prediction 
obtained in the framework of the QED-like effective models can be put under a decisive test against experimental 
data which will be the ultimate means of ascertaining the validity of these scenarios. 



II. EXACT VERSUS LIMITED GAUGE INVARIANCE OF FERMION AMPLITUDES 

In the recent theory of the pseudogap phase of the high-Tc superconductors [4-6], the authors advocated the use of 
the "stringy ansatz" with the inserted Wilson- like exponent of the line integral 

Gr{x-y) =<0\xl;{x)eM-^ J^Mz)dznHy)\0> (2) 

as a viable candidate for the gauge-invariant physical electron propagator (in spite of its being gauge-independent, 
the function Gr explicitly depends on the choice of the contour T which connects the end points x and y). 

Moreover, the authors of Refs. [4-6] argued that the amplitude (2) exhibits the so-called Luttinger-type behavior 

Grix) ^ ^ (3) 

characterized by a positive (albeit varying from one reference to another) value of the anomalous exponent r], provided 

that the contour F is chosen as the straight-line segment between the end points. 

Although the idea of using Eq.(2) for elucidating the behavior of the physical electron propagator has been previously 
entertained, one must be cautioned by the fact that, apart from its relative simplicity, it never received any firm 
justification. A potential problem with both the simple form of the line integral and the particular choice of the 
contour F in Eq.(2) stem from the fact that in the theory of Refs. [4-6] the phase </)(a;) of the (singular) gauge 
transformation '^{x) = 6**^(^^(2^) which relates the physical electron operator '^(x) to the Lagrangian fermion field 
tp{x) is not uniquely defined. 

In what follows, we address some of the misconceptions pertaining to the use of Eq.(2) in the earlier works and 
ascertain the status of the previously obtained results. To this end, we start out by casting Gr{x — y) in the form of 
a functional integral taken over both the gauge and fermion Lagrangian variables 

Gr{x -y) = J D[A]D[i;]D[^Wx)^iy) exp(i ^ A^iz)dz'') exp(i5[V',?, A]) (4) 

Upon integrating the fermions out, the amplitude (4) turns into a functional average over different gauge field config- 
urations 

Gr{x-y) = J D[A]g{x,y\A)exp{iJ^A^{z)dz^')exp{iSeff[A]) (5) 

applied to the inverse Dirac operator g{x,y\A) =< x\[id + A{z) — m]~^\y > computed for a given gauge field 
configuration Af^{z). 

The weight of the average in Eq.(5) is controlled by the efi'ective gauge field action 

Seff[A] = ^ J di^{d^A, - d^A^f + ^Jdxj dyA^{x)Il^Ax - y)A^{y) + • • • (6) 

where the dots stand for the higher order (non- Gaussian) terms produced by fermion polarization (processes of 
"light-light scattering" and alike) which we hereafter neglect, as in all of the previous works on the subject. 

With the transverse fermion polarization Tli^,y{q) = U{q){6i_iv — QnQv/q"^) taken in to account, the gauge field propa- 
gator computed in the covariant A-gauge assumes the form 

q'^+gmiq) q^ 

In order to obtain Eq.(7) one has to complement (6) with a (generally, non-local) gauge fixing term ASeff[A] = 
Jd^JdydM^) < x\[l-g^U{td)/d']\y > d,A,iy)/2g^X. 
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In the case of the conventional (that is weakly coupled, <C 1) QED4, the fermion polarization behaves as 
ll{q) (X q^, and the above gauge fixing term becomes local in the coordinate space. As a result, the ultraviolet (UV) 
divergencies occur at fermion momenta of order the upper momentum cutoff. 

In contrast, QEDj, develops its highly non-trivial behavior at intermediate fermion momenta m <C p ^ A = Ng^ 
controlled by the (this time, dimensionful) coupling constant g [2]. However, adhering to the customary terminology, 
we will refer to this regime as the UV one, for this is where all the logarithmic corrections originate from, whereas 
above the cutoff A no terms oc ln(A/p) can possibly occur. 

In this regime, the gauge propagator (7) is totally dominated by the fermion polarization which, due to the parity 
conserving structure of the reducible four-fermion representation produces no Chern-Simons terms and, in the leading 
1/A'' approximation, reads as 

n(g) = f\/? (8) 

In this way, the number of fermion flavors N becomes the actual parameter controlling perturbation expansion, 
regardless of the strength of the bare coupling g. 

Without making any additional approximations, the function Q{x, y\A) can be expressed in the form of a quantum 
mechanical path integral [11] (see also [12] where this representation was used in a non-perturbative calculation of the 
Dirac fermion propagator in a static spatially random gauge field which is pertinent to such problems as the effect of 
vortex disorder on the quasiparticle properties of rf-wave superconductors or that of dislocations in layered graphite) 



poo i-x{t)=X r-T i n 

g{x,y\A)= dT DxDpe'^°^-' P^ exp{i dr'A^iz) — ) 

Jo Jx(0)=y Jo d,T 



(9) 



where x^{t) is a fermion trajectory parametrized by the proper time t. The spinor structure of the fermion propagator 
is fully accounted for by integrating over the additional variable P(n(r) with the free (matrix- valued) fermion action 

r dx^ 

So[x{t),p{t)] = dr'p^i^ -r- m], (10) 

thus providing a systematic improvement of the celebrated Bloch-Nordsieck model. In the latter, all the spin-related 
effects are ignored which makes this model exactly soluble but restricts its applicability to the infrared (IR) regime 
— m^l <C near the fermion mass shell. In the IR regime, the relevant fermion trajectories contributing to (9) 
deviate only slightly from the straight-path contour F (which coincides with the fermion world line in the case of a 
light cone-like separation between the end points, {x — y)^ = 0). 

One must recognize that the IR regime can only exist if the fermions are massive, while in the case of m = the 
entire region below the upper cutoff A falls into the opposite, UV, regime. In the UV regime, the fermion trajectories 
which dominate the amplitude (9) may strongly depart from the straight-path contour F, for there is no mass term 
to suppress such deviations. 

Combining Eqs.(5) and (9) we represent the amplitude (2) in the form manifesting its gauge invariance 

POO px{t)—x p 

Grix -y)= dT DxDpe'^°^'''P^ < exp(i f d^A^dJl'^'') > (11) 

Jo Jx(0)=y J 

where the brackets stand for the (normalized) functional average over the gauge field with the weight ex.p{iSeff{A]) 
which is to be taken for each trajectory a;^(r) before the sum over all the different trajectories is executed. 

In order to obtain the area integral in the exponent we applied Stokes' theorem to the line integral taken along a 
closed contour which is composed of a trajectory Xh{t) and the straight-line segment F, parametrized as x^^\t') = 
t'{x — y)^/T and traced backwards. 

The averaging procedure generates a sum of all the multi-loop diagrams with no couplings between the fermion 
polarization insertions into the gauge field propagators and the open fermion line which describes the fermion propa- 
gating from y to x. 

Elucidating the structure of the gauge invariant amplitude (2) helps one to appreciate the difference between (11) 
and all the other (^ 7^ 0) members of the family of functions labeled by a continuous parameter ^ (not to be confused 
with the covariant gauge parameter A!) 

_ ^ < Q\i^{x) exp(z(^ - 1) /r A^{z)dznHy)\Q > _ 
^^"^ ^' <0|exp(z^/rA^(z)dz^)|0> 
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h Jm=y <exp(i^/pA^dz'') > 

Amongst others, Eq.(12) also includes the function 

^ „^_ <Omx)^{y)\0> 

""'^^ <0\exp{tJ^A^{z)dz'^)\0> ^ ""^ 

which was invoked by the authors of Ref. [5] who asserted that Gi{x) given by Eq.(13) is identical to the amplitude 
(3) and, therefore, can be used to compute Gr{x) = Go{x). 

The UV anomalous dimension of Eq.(13) can be deduced rather straightforwardly by computing the ratio between 
the wave function renormalization factor determining the anomalous dimension of the ordinary (gauge variant) fermion 
propagator [14] 

< 0\iP{x)i^{y)\0 >(x G{x - y)(A\x - y\)W3^'N){2-3X) (-^4-) 

(hereafter G{x) stands for the free fermion propagator which also yields the bare value of G^{x) for any ^) and the 
Gaussian average of the Wilson line 

< exp(i ^ A^dz^ >= exp[-i ^ < ^ dz!^D^,{zi - Z2)] oc (A|x - y\)(^/-'^K^-^) , (15) 

thus indeed resulting in the anticipated Luttinger-type behavior of Eq.(13) which is characterized by the overall 
positive anomalous dimension 

Notably, this result is free of the gauge parameter A, thus creating the impression that Eq.(13) represents a truly 
gauge invariant function. 

As one immediate objection to the assertion of Ref. [5], one might recall that the very nature of the relationship 
between the electrons and the auxiliary Lagrangian fcrmions implemented via the singular gauge transformation 
< 0|\l'(a;)\l'(j/)|0 >=< 0\il){x) exp[i{^{x) — (f){y)]^/j{y)\0 > implies that the physical electron propagator can only be 
given by a single average over the gauge field, rather than a raiio of such. 

However, except for the case of the original amplitude Go{x), Eq.(12) is represented by a ratio of the two averages 
and, therefore, it can not be truly gauge invariant. Indeed, the gauge fields, over which one averages in both the 
numerator and denominator of Eq.(12), may transform totally independently of one another {Aj^'^ A]/^ + d^f^'"^), 
thus resulting in the overall phase factor exp[i^(/^(a;) — P(x) — f^iy) + P{y))] which only vanishes for / (x) = P{x), 
as long as ^ ^ 0. 

Nevertheless, as shown below, G^{x) turns out to be independent of the gauge parameter A within the class of 
the covariant gauges for any ^. Presumably, it was this limited gauge independence that led the authors of Ref. [5] 
to their conclusion that the (seemingly gauge-invariant) amplitude Gi{x) must be identical to the (truly invariant) 
amplitude Gq{x), for they apparently coincide with each other in the axial gauge {x — y)^An{z) = 0, in which the 
line integral Afj,{z)dz^ vanishes. 

However, should this happen to be valid, the same argument would also apply to an arbitrary G(^{x) which is 
seemingly gauge invariant to absolutely the same extent as Eq.(13) is (that is, provided that one uses the identical 
gauge transformation in both numerator and denominator), and all these functions coincide with each other when 
computed in the axial gauge. 

Nonetheless, one can readily see that G^{x) does bear a non- trivial ^-dependence by simply evaluating the integrand 
in Eq.(12) for certain chosen trajectories. As an example, we consider the case of a square-shaped closed contour 
consisting of four straight segments, each of length |a; — y|. Having all the adjoint sides of the square at right angles 
to each other guarantees one from a possible occurrence of the additional ("cusp") UV singularities [11]. 

Applying the method of dimensional regularization near = 3 we obtain 

< exp(i § d^A^dY.^^ + i^ /r A^dz^) > 
< exp(i^ Jp A^dzi^ > 

exp[-^^^{jp^^r(^)(2 - 0[/i(l - (A - l){d - 2)) - 272(A - l){d - 1)]] 
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where the terms in the exponent proportional to 



' Jo Jo |a-/3|'^-i (3-rf)(2-d) ^ 



and 



' " Jo Jo ip? + /j2)(i+<i)/2 - (3 _ _ \' 1] iiy) 

correspond to the Hne integrals taken along the same and two adjoint sides of the square, respectively (pairs of the 

opposite sides produce no terms oc ln(A|.T; — which is what wc arc after). Despite the fact that the gauge parameter 
A cancels out, as expected, the values of the exponent given by Eq.(17) are markedly different for, say, ^ = and 

In fact, Eq.(17) can cither vanish or grow with increasing \x — depending merely on whether ^ is greater or 
smaller than 2. This observation suggests that for different values of ^ the functions G^{x) are indeed different (it 
would be rather fortuitous if the sum over all the fermion trajectories were ^-independent, despite the difference 
between contributions from each individual trajectory). 

It is worth mentioning that in the case of massive fermions the different functions (12) do become approximately 
identical in the vicinity of the mass shell (|p^ — 4C m^) where the leading functional dependence on the gauge 
field originates from the eikonal phase factor Q{x,y\A) oc cxp(i Jj^A^dz^^). Therefore, the entire IR divergence of the 
ordinary (gauge variant) propagator can be obtained by simply averaging this exponential factor, while any G^(x) 
appears to be totally IR divergence-free. 

This factorization of the exponent of the line integral in the IR regime which had long been known in the case of 
QEDi [13] was extended into the 2>D case in Ref. [9], where it was shown that for any ^ the function (12) possesses 
the same simple pole in the momentum representation (p) oc [p — m)/ [p^ — w?), which corresponds to the universal 
(^-independent) behavior in the coordinate space G^(.t) oc e~'"l'^l. 

In order to avoid confusion we reiterate that the universal IR behavior by no means prevents the different functions 
(12) from having different UV anomalous dimensions, as manifested by Eq.(17). Incidentally, in the massless case it 
turns out to be the UV exponent r/g which controls the power-law decay (3) extending to arbitrarily long distances. 

III. PERTURBATIVE CALCULATION OF FERMION AMPLITUDES 

Having demonstrated the possibility of a non-trivial ^-dependence non-perturbatively, we now present a direct 
perturbative calculation of G^{x). Expanding (12) to second order in Aij,{z) we find three different kinds of correction 

terms 

G^ix -y) = G{x-y)- j dzidzs < G{x - zi)A{z{)G{zx - Z2)A{z2)G{z2 - y) > + (20) 



+(C - l)G{x -y)< A^{zi)dz'^: A^{z2)dz^ > +(1 - / dzi < G{x - zi)A{zi)G{zi -y) A^{z2)dz^ 



> 



The first term in the r.h.s. of (20) corresponds to the first order self-energy correction to the ordinary propagator 
which is given by the leading term of 1/A^ expansion of Eq.(14) 

6,G^{x -y) = -J dz,dz2G{x - zi)rD^Az, - Z2)YG{z2 - y) = G{x - y)^^{'^ - 3A) ln(A|a; - y\) (21) 

The second term in (20) originates from 1/7V expansion of the Wilson line (15) inserted into the numerator and 
denominator of Eq.(12) with the prefactors (^ — 1) and ^, respectively 



62G^{x -y) = {^- \)G{x -y)j^ dz^ ^ dz!^D^,{z, - Z2) = G{x - 2/)^(l - 2^(2 - A) ln(A|x - y\) 



(22) 
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The third term in (20) stems from expanding both G{x,y\A) and the Wilson line in the numerator of Eq.(12) to first 
order in A^{z) 

SaG^ix - J/) = (1 - y dziG{x - zi)rG{zi -y) dz^D^^{zi - z^) « 

« (1 - e) y dzi [G{x - vh'^Gizi -y) + G{x - z^)rG{x - y)] dz^D^^z, - Z2) = 
- OG{x - y)r I dz^, J + 2(1 - A)K ' ^ '^'^^ - 



n^N' ^^^^-^ Jr \l ^\z,-y\^\z,-Z2\^' ' ' ' W- z^? 

= G(x-y)-^{\-i)\\^{h.\x-y\) (23) 

This result was obtained with logarithmic accuracy, and in the course of the calculation we used the coordinate space 
representation of the 3D propagators 



= ^..(-) - ;^2]^[^V + 2(1 - A)^] (24) 



and the following rf-dimensional integrals 



and 



/ 



'^^\x\<^\x-y\'^-^ r(d/2) lyl-^-i ^^^^ 



x<-(x<'-y^)(xi-y'^) _ 2y"^/^T - - j/^J"^ . y"yV . , . 

\x\-i\x - y\-i+^ 3(d-l)r(d/2)^ ^^'^ '\y\'^+^^ ^ ' 



The logarithmic dependence in Eqs. (21-23) results from taking the — > 3 limit of the line integral Jp dz^{z — y)^/\z — 
t/l'^-i = (A|x - y|)3-V(3 -d)^ ln{A\x - y\). 

Combining Eqs. (21-23) together we read off the overall UV anomalous dimension 

^''=35v('^-') ^^'^ 
In particular, for ^ = 1 we reproduce Eq.(16), while for ^ = one obtains 

Thus, despite the fact that the gauge parameter A cancels out, as expected, the functions G^{x) are indeed different 
for different values of the parameter ^. 

At first sight, the non-trivial dependence of Eq.(27) on the parameter ^ (which is linear to first order in 1/N) may 
seem to be in conflict with the fact that all the functions G^ {x) coincide with each other when computed in the axial 
gauge. 

In fact, there would only be a contradiction, if Eq.(12) were truly gauge invariant. Instead, albeit invariant within 
the class of covariant gauges, the function G^{x) may take a different value in the axial gauge, for the latter does not 
belong to this class, for any £, ^ 0. 

On the other hand, the unique property of the function Go{x — y) which is given by a single average instead of 
a ratio of the two does guarantee its true gauge invariance and the possibility to obtain Eq.(28) in arbitrary gauge, 
including both the covariant and the axial ones. In fact, it was shown in Ref. [15] that these two calculations do agree 
with each other, and yet another independent calculation employing the so-called "radial" (Fock-Schwinger) gauge 
{x — yYAf^{x) = [16] provides further support for this conclusion. 

However, the finding that the UV exponent rjof' assumes the negative value suffices to disqualify the "stringy ansatz" 
(2) from being a sound candidate for the physical electron propagator, since, in all of the QED-Vike models discussed 
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so far, repulsive electron-electron interactions arc expected to result in a suppression, rather than enhancement, of 
any amplitude describing propagation of physical electrons. 

The above analysis can be easily extended to the case of the conventional weak coupling QED4 which demonstrate 
that the situation in 3D is not at all exceptional. Instead of Eqs. (21-23) we now get 

d,G^{x -y) = -G{x - y)^X\n{K\x - y\) (29) 

52G^{x -y)= G{x - - 2^(3 - A) ln(A|:r - y\) (30) 

5^G^{x -y) = G{x - y)£^{l - OAln(A|a; - y\) (31) 
and, instead of Eqs. (24-26) , we use the formulas 



|a.|d|a._y|d-2 r(d/2) |y|<^-2 



(33) 



/ 



x^{x^ -yP){x^ -y^) _ TT'^/^ ^ y^S^-* - - y^5"^ m_ 9^^!^1 

\x\'i\x-y\<i ~ 2{d-2)T{d/2Y |y|''-2 + ^> \y\<i ^ ^"^^^ 



Combining Eqs. (29-31) together we obtain the UV anomalous dimension of G^{x) in the 4D case 

^r = ^(2^-l) (35) 

Interestingly enough, for ^ = and ^ = 1 the values of Eq.(35) differ only in their sign, and the anomalous dimension 
Vio^ = —Sg^/SiT^ is again negative. 

IV. ALTERNATE FORMS OF PHYSICAL ELECTRON PROPAGATOR 

The apparent problem with the unphysical decay of the conjectured form of the physical electron propagator (2) 
which appears to be slower than in the case of non-interacting electrons compels one to explore alternate proposals. 
To this end, one may be able to benefit from the previous studies of the conventional QED^ where the problem of 
constructing gauge invariant asymptotical quantum states with quantum numbers of an electron (" dressed charges" ) 
was long recognized, and attempts to explicitly solve it have long been under way. 

A (non-local) composite operator creating such a state can be sought in the form 

^{x) = exp[i J dyrix - y)A^{y)]i^{x) (36) 

where the vector function f^{x) is subject to the condition d/if^^x) = 5{x) imposed by the requirement of gauge 
invariance. 

In the original proposal which was put forward by Dirac in his pioneering work of the fifties, the physical electron 
was identified with the operator 

*^(a:) = eMi^mx) (37) 

Indeed, the expectation value < 0|4'£)(a;)A^(y)4'|j(a;)|0 >= 5^o/47r|£ — y\ describes the Coulombic electrostatic field, 
thus suggesting that Eq.(37) provides a plausible representation of a static charge. 
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Moreover, it was found that the Fourier transform of the gauge-invariant amplitude < 0\'^ d{x)^ D{y)\0 > is IR 
finite at the single point = {m, 0) on the mass shell [17]. 

It was also argued in Ref. [17] that, once created, the coherent state 4'|5(x)|0 > remains stable, unlike that produced 
by the operator "^rix) = exp(i Ao{t')dt')ip{x) whose propagator < 0\'^Tix)^T{y)\0 > is given by Eq.(2) where 
{x - y)^ oc (5^0- _ 

Furthermore, the amplitude < d{x)^ D{y)\Q > was found to undergo multiplicative UV renormalization, fea- 
turing a positive UV exponent 

Vjf = g^Sn' (38) 
The authors of Ref. [17] were also able to generalize the construction (37) onto the case of a moving charge 

*'(^* = '"■x-lpOT^^"'^'***^* '''' 

where = (1,0) is the time-like unit vector and t;^ = {0,v) is the velocity of the charge. The formula (39) was 
obtained by requiring that, in accord with the anticipated physical interpretation of Eq.(39), the electromagnetic 
field < 0\'^ v{x)Af^{y)'i'l{x)\0 > associated with the coherent state created by this composite operator reproduces the 
classical Lienard-Wiechert potentials. 

It was also found that the corresponding propagator 

G.ix - y) =< Omx) exp[^ _ ^ ^^^^J d,A.]i,{y)\0 > (40) 



is both IR-finite at = m{l, v)/\/l — v'^ and multiplicatively UV-renormalizable, its anomalous dimension being 
explicitly u-dependent [17] 

The combination of all the befomentioned properties makes Eq.(39) the best (up-to-date) available candidate for the 

physical electron operator in the conventional (massive) QED/^. 

It is worth mentioning that Eq.(41) coincides with the UV dimension of the 4Z) "stringy" amplitude (2) only in the 
unphysical limit t; — > oo, while its u = value yields Eq.(38). 

By construction, Eq.(41) can not be applied directly to the massless, let alone the 3D massless, case which corre- 
sponds to the limit v ^ 1 where the UV dimension (41) increases towards infinitely high positive values. 

Albeit not being immediately applicable to the case of m = 0, the calculation of Gy{x) carried out in Ref. [17] 
suggests that the logarithmic growth of (41) gets effectively cut off at ma.T(|l — ?;|,1/A|.t|). Therefore, it is not 
inconceivable that the massless counterpart of Eq.(40) may exhibit a faster than a power-law decay 

Gphys (x) (X exp ( — mnst In^ ( A | x | ) ) , (42 ) 

where the constant is proportional to either 1/N (3D) or (4D), thus placing the effective QED-Mkc theories into 
the class of " super-Luttinger" models, alongside the one-dimensional metals with unscreened Coulombic interactions 
where Gphys{x) oc exp(— const ln^/^(A|a;|)) [18]. 

In light of this possibility, the previous claims of discovering the Luttinger-like behavior in the framework of the 
Q-E-Da-theory of the pseudogap phase [4-6] may need to be prepared to handle a potentially much stronger suppression 
of the physical amplitudes in order to reconcile such a behavior with the available photoemission, tunneling, and other 
experimental data. 



V. DISCUSSION 



Evidently, one could escape the whole problem with the singular massless limit, if the Dirac fermions acquired a 
finite mass via the mechanism of chiral symmetry breaking [19]. In the context of the QED^^ theory of the pseudogap 
phase, such an intrinsic instability is expected to result in the onset of a spin (SDW) and/or charge (CDW) density 
wave [7,20], provided that such a transition does occur for the physical number of fermion flavors [N = 2). 
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Thus far, there have been only prehminary conclusions drawn on the role of the strong spatial anisotropy of the 
quasiparticle dispersion (e = \/vfqf + t^j^l with Vi/v2 > 10) in the high-Tg cuprates. 

Although, in accord with general physical expectations, this anisotropy was found to decrease upon IR renormal- 
ization, [20,21], there has been no consensus reached, regarding universality (or a lack thereof) of the critical number 
of flavors Ncr, much less regarding the possibility that the spatial anisotropy can drive this critical number to the 
values in excess of 2. 

To this end, it may be helpful to complement the renormalization group approach adopted in Refs. [20,21] by a direct 
analysis of the strong coupling regime where the gauge field propagator (7) is dominated by the fermion polarization 
n(g). In the case of anisotropic fermion dispersion, the latter becomes a function of the spatial components of 
the momentum q multiplied by the corresponding velocities ViQi, thus allowing one to scale such factors out in the 
Schwinger-Dyson equation for the fermion mass [19] and suggesting that N^r may even remain the same as in the 
isotropic case. 

In any event, the effect of the spatial anisotropy of the fermion dispersion may be rather different from that of the 
Lorentz non-invariant (e.g., pure Coulombic) couplings which can cause N^r to decrease [9,10], as compared to the 
value obtained in the Lorentz-invarint case. 

One should, however, be alerted by the fact that even in the perfectly Lorentz-invariant case the status of the 
analyses based on the Schwinger-Dyson equation (which, in this particular case, suggests N^r > 3) remains somewhat 
unclear. In this regard, it was pointed out in Ref. [24] that the Schwinger-Dyson equation may systematically 
overestimate the actual critical number of flavors which might, in fact, be as low as 3/2. 

It is, therefore, quite remarkable that, as an alternate route, the intrinsic SDW instability of the pseudogap phase 
of the cuprates can instead manifest itself through a divergent staggered spin susceptibility [22,23] 

XQ{q) oc (y^)i-64/(3^^iV) (43) 

where the momentum q corresponds to the deviation from the SDW ordering vector Q = (it, tt). 

In order for this behavior to occur, the presence of massless neutral fermions is rather necessary than problematic. 
Curiously enough, the susceptibility (43) displays the negative anomalous dimension which is exactly twice the value 
(28). 

Although this observation does not necessarily imply that the staggered spin susceptibility given by the aver- 
age xq{^ ^ v) =< Q{x,y\A)Q{y,x\A) > is directly related to the product of the two averages \Gn{x — x< 
Q{x,y\A)ex.p{—i J^A^dz"^) >< G{y,x\A) exp{i J^Af^dzf^) >, it may help one to elucidate the real physical meaning 
of the negative anomalous dimensions exhibited by the gauge invariant amplitude (2). 

It is possible that, albeit being unfit for the role of the physical electron propagator, the amplitude (2) may still 
bear some important information about the vertex functions which also determine the behavior of the gauge-invariant 
susceptibilities. 

As far as the experimental status of the QED^ theory of the cuprates is concerned, it is presently unclear whether 
it can accomodate the phenomenon of time reversal symmetry breaking, as suggested by the recent experiment [25], 
without generating the Chern-Simons term in the fermion polarization. 

In contrast, broken time reversal symmetry seems to be an intrinsic feature of some alternate approaches to the 
pseudogap phase which, instead of the phase fluctuations of the local parent 2 2 -wave order parameter [4-7], 
focus on other bosonic collective modes, such as emergent d-symmetrical CDW [26] or incipient secondary pairing 
{dx:2-y2 ^ dx2-y2 +idxy) [27,28]. 

To this end, it is worth mentioning that the recent experimental reports of a genuine quantum-critical behavior 
in the Ca-doped cuprates [29] are in quantitative agreement with the predictions [28] based on the theory of Ref. 
[27]. Notably, in this alternate (Higgs- Yukawa- type) theory of the cuprates the anomalous dimension of the electron 
propagator was found to be positive [28,30]. 

We conclude by stressing that, regardless of the status of the QED3 theory of the cuprates itself, the problem of 
constructing the gauge invariant fermion propagator in the eff'ective massless QED-like theories, thus far, has received 
a lesser attention than it, arguably, deserves. 

This problem needs to be settled before one can start drawing solid, rather than wishful, conclusions about the true 
behavior in these as well as other gauge field models, including non-abelian and discrete symmetry (say, Zn) ones. 

In the present paper, we undertook an attempt to clarify this issue which may have already resulted in a widespread 
confusion, in particular, as far as the QED^ theory of the pseudogap phase of the cuprates is concerned. We 
demonstrated that the previously proposed ansatz ehxibits a clearly unphysical behavior and, therefore, needs to be 
modified. In the course of this analysis we established an interesting property of "partial gauge invariance" of the 
family of functions (12) with ^ ^ which, we believe, might have been partly responsible for the erroneous conclusions 
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drawn in [5]. We also conjectured an alternate form of the physical electron propagator which decays with distance 
faster than any power-law, thus challenging the interpretation of the experimental data based on the idea of the 
Luttinger-like behavior in the cuprates proposed in Refs. [4-6]. 
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